Probability distribution of the order parameter 
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The probability distribution of the order parameter is exploited in order to obtain the criticality of 
magnetic systems. Monte Carlo simulations have been employed by using single spin flip Metropolis 
algorithm aided by finite-size scaling and histogram reweighting techniques. A method is proposed to 
obtain this probability distribution even when the transition temperature of the model is unknown. 
A test is performed on the two-dimensional spin- 1/2 and spin-1 Ising model and the results show 
that the present procedure can be quite efficient and accurate to describe the criticality of the 
system. 



o 

B 



-a 
c 

o 
o 



CO 
CN 

o 

o 

-i— > 

c3 



i 

C 

O 
o 



X 



PACS numbers: 05.10.Ln, 02.50.Ng, 64.60.-i 



I. INTRODUCTION 

The order parameter distribution function has been 
proved to be apowerful tool for studying not only mag- 
netic systems yj, y, Li La LI , but also the liquid-gas critical 
point 0| , the critical point in the unified theory of weak 
and electromagnetic interactions 0, and the critical point 
in quantum chromodynamics 8J . For the specific case 
of magnetic systems the order parameter can be chosen 
as the magnetization per spin, namely m — jj 2f=i 
where N is the total number of spins and Si is the spin 
at site i. In finite-size systems, the magnetization m is 
a fluctuating quantity, characterized by the probability 
distribution P(m) 0)0- In Ising- like models undergo- 
ing a second-order phase transition it is known that at 
temperatures lower than the critical temperature T c , the 
distribution P(m) has a double peak, centered at the 
spontaneous magnetizations +m and — m. At temper- 
atures greater than T c , P(m) has a single peak at zero 
magnetization, and exactly at T c a double-peak shape is 
observed^. Analogously to the usual finite-size scaling 
assumptions |9j, one then expects that, for a large finite 
system of linear dimension L at the critical point, P(m) 
takes the form 



P(m) = bP*(rh), 



(1) 



where b = b^L^I 1 ' ', [3 and v are critical exponents, 
77i = bm, bg is a constant, and P*(rh) is a universal scal- 
ing function, normalized to unit norm and unit variance. 
Scaling functions, such as P* , are characteristic of the 
corresponding universality class. Systems belonging to 
the same universality class share the same scaling func- 
tions. Thus, from the precise knowledge of P*{rh) one 
can characterize critical points and also identify univer- 
sality classes. This is what has been done so far in the 
literature, with the distribution for the spin-1/2 Ising 



model being the standard P* functionlj, G] for this uni- 
versality class. For instance, it is shown in Fig. ^ the 
normalized distribution P*(m) for the two-dimensional 
spin-1/2, spin-1, and spin-3/2 Ising model at critical- 
ity. Simulations have been done on square lattices with 
L = 32 at the exact T c for spin-1/2, at T c = 1.6935 for 
spin-1, according series expansions 101 and Monte Carlo 
simulations and at T c — 3.28794 for the spin-3/2 
model ■ The universal aspect of these systems can be 
easily noted. 




FIG. 1: Scaling function P*(fh) for the two-dimensional spin- 
1/2, spin-1, and spin-3/2 Ising model on square lattices with 
L = 32. Simulations were performed at the exact T c for spin- 
1/2, at T c = 1.6935 for spin-1, according series expansions 
[iq and Monte Carlo simulations £3, and at T c = 3.28794 
for the spin- 3/2 model HJ. The error bars are smaller than 
the symbol sizes. After Ref. [Tl| 

Monte Carlo simulations seem to be the most effec- 
tive method to obtain results as those shown in Fig^ 
where the probability distribution P(m) corresponds to 
the fraction of the total number of realizations in which 
the system magnetization is m, i.e., 
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where N m is the number of times that magnetization m 
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appears and Nmcs is the total number of Monte Carlo 
steps. To compute the normalized distribution P*(m) via 
Eq. one has to evaluate the pre-factor b. This can 
be easily done by noting that b = l/cr, where a is the 
square root of the magnetization variance (a 2 — (m 2 ) — 
(m) 2 ). Thus, one obtains the universal function P*(rh) 
by simply rescaling the magnetization and by using Eq. 

09- 

In general, the probability distribution is used for 
studying models in which the critical temperature or even 
the distribution function is exactly (or high-precisely) 
known. That is in fact what has been done in the study 
of several systems. When this distribution, as well as 
the critical temperature and critical exponents, are not 
known, one can of course do first a canonical simulation 
in order to get the critical values (universal and non uni- 
versal) and compute, afterwards, the desired distribution. 
The present approach is different from this conventional 
one in the sense that it does use the order parameter 
distribution itself in order to obtain the criticality of the 
system. The procedure, as well as the results obtained 
for the spin-1/2 and spin-1 Ising model, are discussed in 
section II and the conclusions are presented in the final 
section. 



II. APPROACH AND RESULTS 

We have performed extensive Monte Carlo simulations 
(up to 10 7 - 10 s Monte Carlo steps per spin after 2.0 — 
5.0 x 10 4 steps for thermalization) on square LxL lattices 
with periodic boundary conditions for systems of length 
12 < L < 64. For a given L, the simulation ran at a fixed 
temperature, evolving according the standard Me trop olis 
algorithm. A histogram reweighting technique [12l Il3| 
was used to obtain thermodynamic information in the 
vicinity of the simulated temperature. 

Let us first discuss the spin-1/2 Ising model. Figure 
121 shows the distribution P* as a function of the nor- 
malized magnetization m for temperatures different from 
the critical value T c . As expected, one can see that for 
a temperature lower than T c (Fig. the maximum 

value of P* increases when the lattice size L increases, 
while for a temperature greater than T c , an increase of 
L leads to a decrease of the corresponding peaks of P* 
(see Fig. |2t>). In other words, suppose we have a distri- 
bution function P*(rh) for a given L (say for example, 
L = 16) at a fixed temperature Tl = i6- If Tl=16 < T e , a 
similar distribution will be obtained for a bigger lattice 
(e.g., L = 64) at a different temperature Tl=64 such that 
T L= i 6 < T L=64 < T c . Analogously, if T L=16 > T c , we will 
have 2z, = i6 > Tl = 64 > T c . This fact suggests a mech- 
anism to obtain the critical temperature, as well as the 
exponent v and the universal distribution, for the system 
under study. Table [fl shows the temperatures of several 
lattice sizes we have used to obtain the distributions dis- 
played in Fig. |3J These temperatures were evaluated 
as follows. For L = 64 and a given temperature, for 
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FIG. 2: Scaling function P*(rh) for the spin-1/2 Ising model 
with L — 16, 32, and 64 at a fixed temperature T: (a) lower 
than T c (T = 2.2472) and (b) greater than T c (T = 2.2831). 
The error bars have been ommited for clarity. 



TABLE I: Temperature for different lattice sizes at which the 
distribution P*(m) for L = 16 — 48 is the same as that ob- 
tained for L = 64 at the shown temperatures (spin-1/2). Er- 
ror in parentheses affects the last digits. The second and third 
columns correspond to temperatures greater than the critical 
one, and the two following columns correspond to tempera- 
tures lower than the critical one. The last column represents 
the data when P*(rh) for L — 64 is obatined at T c . 

Size Temperature (in units of J/fcs) 

16 2.3923(11) 2.3272(8) 2.2477(8) 2.1901(10) 2.27221(52) 

20 2.3666(8) 2.3154(8) 2.2502(8) 2.2036(10) 2.27092(52) 

24 2.3502(8) 2.3073(5) 2.2528(5) 2.2134(7) 2.27015(52) 

32 2.3288(5) 2.2973(5) 2.2563(5) 2.2262(5) 2.26963(26) 

48 2.3089(5) 2.2878(5) 2.2604(5) 2.2399(5) 2.26937(26) 

64 2.2989 2.2831 2.2624 2.2472 2.269184 



instance T = 2.2989 in Table DJ we compute the corre- 
sponding probability distribution of the order parameter, 
which will be the "reference" distribution. For other val- 
ues of L, we search for the temperature Tl which gives a 
distribution equivalent to the reference one. In this way, 
we obtain the data shown in the second column of Table 
[J] Taking a different reference distribution, obtained at a 
different temperature for L — 64, we have another set of 
Tl, and so on. All the distributions so obtained are de- 
picted in Fig. |3 It means that each curve in Fig. |3|is in 
fact a superposition of six different distributions taking 
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at the temperatures given in Tabled 

2.5 r 




FIG. 3: Normalized distribution P*(rh) for systems with lat- 
tice sizes and temperatures shown in Table [I] Each curve is 
a supperposition of six different distributions taking with the 
data from this table. 

Since one expects that the difference \Tl — T c \ scales as 
L^ 1 /" , where u is the correlation length critical exponent, 
a finite-size scaling analysis can be done to estimate the 
critical values of the infinite system. In Fig. we have 
a plot of Tl vs. L -1 /", with v = 1, using the values of 
the first five columns of Table which confirms the exact 
exponent v = 1 and gives T c — 2.267(2). Another choice 
for this model is, of course, the corresponding critical 
temperature. In this case (very close to T c , last column 
of Table QJ, however, it is known that \Tl — T c \ scales as 
£,-(1+0)/^ where 8 is the correction to scaling exponent 
0. In Fig. 0Jd, we plot the estimates Tl as a function of 
£-(1+0)/" w ith v = 1 and 8 — 2 14]. Linear regression 
gives T c — 2.2693(1) for the infinite system, which is in 
fact quite close to the exact one. 

In order to measure the applicability of the present 
mechanism for obtaining the transition temperature from 
non-universal distributions, we also study the spin-1 Ising 
model. To have an idea of the value of T c , one just per- 
forms short simulations in a range of temperatures to 
check whether the probability distribution P(m) has sin- 
gle or double peak. Then, one proceeds according to the 
same manner already discussed for the spin-1/2 case. We 
fix the temperature and verify how the peaks of the dis- 
tribution change if the lattice size L increases. Figure 
shows the distributions obtained for lattice sizes L = 12, 
16, 24, and 32 at two different temperatures: T = 1.660 
and T — 1.720. In the former case an increasing lattice 
size leads to increasing peaks, while in the latter case the 
height of the peaks decreases when the lattice becomes 
larger. Thus, we conclude that the transition tempera- 
ture is between 1.660 and 1.720, and hence we perform 
longer simulations in this temperature range. 

The procedure now is the same as that we have done 
for the spin-1/2 model. We fix the temperature and com- 
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FIG. 4: (a) Temperature Tl as a function of L -1 ^" with v = 1. 
Different symbols correspond to different choices of the ref- 
erence distribution (from top to bottom we have the results 
from the second to fifth columns of Table [TJ. Error bars are 
smaller than the symbol sizes, (b) Temperature Tl as a func- 
tion of L" (1+e)/ " with v = 1 and 9 = 2 taking the data of the 
last column of Table Q 



pute P*(rh) for the lattice with L = 32 (reference distri- 
bution). For a different L, we search for the temperature 
that gives a distribution equal to the reference one. Ta- 
ble [H] shows the temperatures so obtained. For each set 
of temperatures (which corresponds to each column of 
Table [nj, we plot Tl vs. L" 1 ^ and vary the exponent 
v until we get a straight line. Thus, each column of Ta- 
ble [H] gives an independent estimate of v and also of the 
critical temperature T c . Figure illustrates this proce- 
dure. By taking the mean value of these quantities, one 
obtains v = 1.0(1) and T c = 1.6933(16), where the latter 
agrees well with the value T c = 1.6935(10) 10]. After we 
have evaluated T c , we ran a longer simulation on a larger 
lattice to determine, by this way, the universal distri- 
bution P*(rh). Figure [7| shows the distribution P*(rh) 
on a L — 64 lattice for spin-1/2 and spin-1 models, and 
confirms the fact that both systems belong to the same 
universality class, as already expected. 



III. CONCLUSIONS 

The present approach, using just the order parame- 
ter distribution, seems to be a robust way to obtain the 
criticality of magnetic systems, regarding its universal 
and non-universal aspects. There is a clear distinction 
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FIG. 5: Scaling function P*(rh) for the spin-1 model on square 
lattices with L = 12, 16, 24, and 32 at a fixed temperature 
T: (a) lower than T c (T = 1.660) and (b) greater than T c 
(T = 1.720). 



TABLE II: Temperature for different lattice sizes at which the 
distribution P*(m) for L — 12 — 24 is the same as that ob- 
tained for L = 32 at the shown temperatures (spin-1). Error 
in parentheses affects the last digits. The second and third 
columns correspond to temperatures greater than the critical 
one, and the two following columns correspond to tempera- 
tures lower than the critical one. 

Size Temperature (in units of J/ks) 
12 1.598(4) 1.652(2) 1.708(2) 1.765(1) 
16 1.624(2) 1.664(2) 1.705(2) 1.747(1) 
24 1.649(2) 1.675(1) 1.702(1) 1.730(1) 
32 1.660(2) 1.680(1) 1.700(1) 1.720(1) 



between the finite-size behavior of P* close to the criti- 
cal temperature (scaling with or away from it 
(scaling with L -1 /"), as depicted in Figs. ^and0 It also 
seems, at first sight, that there is a limitation regarding 
the size of the lattices which could be considered. For in- 
stance, we have used here lattice sizes which are smaller 
than that from the reference distribution. Nevertheless, 
this limitation is not so drastic since in the spin-1 model 
we considered the reference distribution for L = 32 and 
with lattices smaller than this value the results prove 
quite accurate. One can, of course, consider lattices 



FIG. 6: Temperature Tl as a function of L' 1 ^ (spin-1). Dif- 
ferent symbols correspond to different choices of the reference 
distribution (see Table ITTI . The values of v that give the best 
linear fit were (from top to bottom): v — 0.9(1), v — 0.9(1), 
v — 1.1(1), and v = 1.1(1). Error bars are smaller than the 
symbol sizes. 
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FIG. 7: Universal function P*(m) for lattice size L = 64 at 
temperatures obtained in this work: T = 2.2693 for spin-1/2 
and T = 1.6933 for spin-1. 



larger than that of the reference distribution. We feel, 
however, that reweighting the distribution is more eas- 
ily done for smaller systems. Application of the present 
procedure to other models (pure and random), as well as 
to multicritical behavior, will be very welcome; some are 
now in progress. 
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